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Chapter 1 


Lotka- Volterra Dynamical Systems 


Stephen Baigent 


Department of Mathematics, UCL, Gower Street, London WC1E 6BT 
steve. baigent@ucl.ac. uk 


Lotka-Volterra systems are used to introduce in a simple setting a num- 
ber of dynamical systems techniques. Concepts such as omega limit 
sets, simple attractors, Lyapunov functions are explained in the context 
of Lotka-Volterra systems. We discuss LaSalle’s Invariance principle. 
Monotone systems theory is also introduced in the context of the Lotka- 
Volterra systems. 


1. Introduction and scope 


The Lotka-Volterra equations are an important model that has been widely 
used by theoretical ecologists to study the implications of various inter- 
actions between members of a population in a fixed habitat containing a 
number of distinct interacting species. They are by no means the most real- 
istic of such ecological models, but they are arguably the simplest since the 
highest order terms they involve are quadratic, and therefore they feature 
the next level of complexity up from linear. As we shall see, even amongst 
differential equations with quadratic terms, they have a very special form 
which makes them amenable to well-known mathematical techniques from 
standard linear algebra, convex analysis, and dynamical systems theory. 

To set the scene, we write the Lotka-Volterra equations in the revealing 
form: 

t; n 
7 

aba (1) 
where n is the number of distinct species, x; is the population density of 
the ith species, and r;,a;; are all real numbers, possibly zero, and here 
assumed to be independent of time. Multiplying each equation through 
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by x; shows that indeed the equations are quadratic, but when written as 
above we see that the net population growth per individual per unit time 


(“;/a2;) is linear in the population densities x = (a1,...,2n). 
A good starting point in the study of the dynamics of (1) is to first 
locate steady states; that is, points x* where t; = 0 for each i = 1,...,n. 


Of special interest, since they model one scenario where all species can 
coexist, are the so-called interior steady states. These satisfy x7 > 0 for 


each i = 1,...,n and so are obtained by solving the linear system 
n 
ri+ Y ayti, i=1,...;n. (2) 
j=1 


Recall that the r;,a;; may be of any sign or zero. As we shall see, deter- 
mining when (2) has a unique solution x} > 0 for each i = 1,...,n relies 
heavily on linear algebraic techniques. All other steady states involve at 
least one density vanishing; that is at least one species is extinct. Such 
steady states are determined by investigating the linear system (2) with all 
possible subsets of {a;}"_, set to zero. 

The main virtue of model (1) is that it enables us to study on paper 
or on the computer the outcome of any set of interactions between the n 
species, and they are the simplest model that enables us to do so. The 
type of interactions we refer to are split into two categories: intraspecific 
(the effect of one member of a species on another member of the same 
species), and interspecific (the effect of a member of one species on a mem- 
ber of another species). The strength of the interactions are encoded in 
the parameters aij, which are usually assembled into the n x n interaction 
matrix A = ((a;;)). The parameters r; determines the intrinsic growth 
rate per individual of species i which would be observed if intraspecific and 
interspecific interactions were absent. Here we will not be concerned with 
the precise details of the ecological or environmental mechanisms that con- 
tribute to the value of each of the parameters, as we are more interested in 
the effects of the signs and magnitudes of each parameter on the qualitative 
behaviour of (1). 

We shall however, link the signs of parameters to types of interactions. 
For example, r; > 0 says that the environment intrinsically favours the 
growth of species 7, whereas r; < 0 signals a risk of extinction for that 
species unless the presence of another species promotes its growth. An 
example of the latter case is where a predator will go extinct in the absence 
of its prey and a suitable substitute food source. 

Much of the material will apply to the most general form of the Lotka- 
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Volterra model (1). Existence of interior steady states will be investigated 
and their local stability studied. 


2. Lyapunov methods for Lotka-Volterra Systems 


2.1. Some basic dynamical systems results 


Sime notation first: R>ọ = {x > 0}, R>o = {x > 0}. We will always 
assume that parameters are such that the differential equations (1) generate 
a semiflow y; : R&, > R&,: Va € R&, and s,t > 0, 


(1) Yo(x) = z; 
(2) vt(ps(x)) = Yt45(2); 


Let U C R” be open. 


Definition 1 (Orbit). The (forward) orbit of x € U is the set OT (x) = 
{yi (x) :t > O}. 


Definition 2 (Steady state). A steady state of t = f(x) isa point x € U 
for which f(x) = 0. 


Definition 3 (Forward invariant set). A set S C U is a forward invari- 
ant set for p+ if whenever x € S we have y(x) € S for all t > 0. 


Definition 4 (Invariant set). When ys is a flow (i.e. also defined for 
t < 0), the set S C U is an invariant set for y, if whenever x € S we have 
plx) € S for allt € R. 


One important use of invariant sets is captured by the following result:? 


Theorem 1. Let S C R” be homeomorphic to the closed unit ball and 
forward invariant for the flow of t = f(x). Then the flow has a steady 
state x* € S. 


Hence one way of showing the existence of at least one steady state in a 
compact simply-connected subset of R” is to show that all orbits enter that 
set (so that it is forward invariant). 

The Heine-Borel theorem states that a subset of R” is compact if and 
only if it is closed and bounded. The key tool for studying the convergence 
of orbits is the Omega limit set. This is the totality of all limit points of the 
forward orbit through a given point. To prove that an orbit is convergent 
to a steady state, one needs to show that its omega limit set consists of 
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a single point, namely that steady state. Other interesting limit sets are 
attracting limit cycles, periodic orbits, attractors, etc. 


Definition 5 (Omega limit point). A point p € U is an omega limit 
point of x € U if there are points pr (x), pta (x),... on the orbit of x such 
that tk > co and yz, (x) > pas k > œœ. 


Definition 6 (Omega limit set). The omega limit set w(x) of a point 
x E€ U under the flow y, is the set of all omega limit points of x. 


There is a similar construct for when y; is defined backwards in time, such 
as when it is a flow: 


Definition 7 (Alpha limit point). A point p is an a limit point for the 
point x € U if there are points ps (£), y1,(#),... on the orbit of x such that 
tk > —oo and p, (x) > pas k > oo. 


Definition 8 (Alpha limit set). The alpha limit set a(x) of a point x € 
U under the flow vy; is the set of all alpha limit points of x. 


Lemma 1 (Properties of Omega limit sets). 


(1) w(x) is a closed set (but it might be empty); 
(2) If O+(a) is compact, then w(x) is non-empty and connected; 


(3) w(x) is an invariant set for pr; 
(A) fy © O* (i) them wly) = w(x); 
(5) w(x) can be written as 


w(x) =(){es(2):s> th = () Ore), 


t>0 t>0 


where A is the closure of A. 
For a proof see, for example, reference 4. 


Example 1. z = 1 has the flow y(x) = x + t. Given any x € R and any 
sequence tg — 00, Y, (£) oo and hence w(x) is empty. On the other 
hand, for « = az the flow is y(x) = ez, so that y;, (x) = e“*x > 0 as 
tk > œ if a < 0 giving w(x) = {0} and clearly y,(0) = 0 so w(x) is indeed 
invariant. But if a > 0 the set w(x) is empty. 


As another example, take 
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Example 2. 


t =g- y- grle? +y’) 
ý =z +y- ylz? +y’). 8) 


By multiplying the first equation by x and the second by y and adding we 
obtain, after setting r = \/x2 + y? and simplifying, * = r — r3. The set 
r=lie. S= {(x,y) : x? +y? = 1} is an invariant set and (x,y) = (0,0) 
is the unique steady state. It is not difficult to see that any orbit is either 
the unique steady state (0,0), the unit circle, or a spiral that tends towards 
the unit circle. If (x,y) # (0,0), w((x,y)) = S, and otherwise w((0,0)) = 


{(0, 0)}. 
Problem 1. Find the omega limit sets for the predator-prey system on Ro 


t=x(1— x+y) 
ý =y(-1 -y+ zx). 


The many practical uses of the omega limit set is typified by the fol- 
lowing result. Note that + = 1/x with z(0) > 0 satisfies t —> 0 as t > ~w, 
but the unique forward orbit x(t) = \/2t + 7(0)? > œo as t > oo does not 
converge to a steady state. However, we do have: 


Lemma 2. Suppose that f : R” — R” is continuously differentiable with 
isolated zeros. If x : Rso > R” is a bounded forward orbit of t = f(x) 
such that z(t) + 0 as t > œ, then x(t) > p for some p as t > co where 
f(p) = 0, i.e. x converges to a steady state. 


Proof. Let the orbit pass through zo. O+(xo) is bounded and hence com- 
pact, so w(zo) is compact, connected and nonempty. For p E€ w(xo) there 
exists a sequence t — co as k — oo such that x(t.) > p as k > oo. 
By continuity 0 = limpoo (tk) = limgoo f(a(te)) = f(p), so that p is a 
steady state. Thus w(zo) consists entirely of steady states. Since w(xo) is 
connected, and the steady states are isolated, w(xo) = {p}. 


2.2. Stability 


Definition 9 (Lyapunov stability). A steady state x* is said to be Lya- 
punov stable if for any e€ > 0 (arbitrarily small) 3 6 > 0 such that Vio with 
|a* — zo| < 6 we have |y(xo,t) — x*| < e for all t > 0. 


A steady state is said to be unstable if it is not (Lyapunov) stable. 
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Definition 10 (Asymptotic stability). A steady state x* is said to be 
asymptotically stable if it is Lyapunov stable and 3 p > 0 such that Vo 
with |xz* — zo| < p we have |p(xo, t) — 2*| > 0 as t > co. 


For example, in the system ¢ = ~g — y + x(x? +y?), ý = x — y +ylr? +y’), 
the origin is locally asymptotically stable (we get * = —r+r° by using polar 
coordinates). For a simple harmonic oscillator in the form of a pendulum, 
the pendulum resting vertically downwards is Lyapunov stable but not 
asymptotically stable unless there is damping such as air resistance. The 
upward vertical state of the pendulum is an example of an unstable steady 
state. 


Definition 11 (Basin of attraction). The basin of attraction B(x*) of 
a steady state x* € U is the set of points y € U such that p(y) > x* as 
t> oo. 


Definition 12 (Global stability). If B(x*) = U then 2* is said to be 
globally asymptotically stable on U. 


Problem 2. Consider the logistic equation t = x(1 — x). Find all forward 
invariant and invariant subsets of R>o and obtain the basin of attraction 
of the positive steady state. 


3. Ecological Systems 


Consider the model 


peni ea (4) 
where each f; : R” — R” is C!. Suppose that z(0) = (xo1,---,; on) has 
zok = 0 for k € J C {1,...,n}, so that some species are initially absent. 


Then these species are absent for all time for which the solutions exist: 


Theorem 2. For the model (4) the coordinate azes and the subspaces 
spanned by them, and R&,, are all forward invariant. 


In other words populations that start nonnegative remain nonnegative. 
Populations starting positive cannot go to zero in finite time. 


4. LaSalle’s Invariance Principle 


We start with a basic result for Lyapunov functions (e.g. page 127 in 
reference 13): 
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Theorem 3. Let U C R” be open and f : U > R be continuously differ- 
entiable and such that f(xo) = 0 for some zo E€ U. Suppose further that 
there is a real-valued function V : U > R that satisfies (i) V (xo) = 0, (ii) 
V(x) > 0 for x € U \ {zo}. Then if (a) V(x) := VV (x) - f(x) < 0 for all 
x € U then zo is Lyapunov stable; if V(x) < 0 for all U \ {x9} then zo 
is asymptotically stable on U; (c) if V(x) > 0 for all x € U \ {xo}, £o is 
unstable. 


This is a powerful theorem, but there is a useful generalisation of it which 
caters for when V~1(0) is not an isolated point. 


Theorem 4 (LaSalle’s Invariance Principle). Let « = f(x) define a 
flow on a set U C R”, where f is continuously differentiable. Suppose 
V : U —> R is a continuously differentiable function. Let Q be the largest 
invariant subset of U. If for some bounded solution x(t,29) with ini- 
tial condition x(0,29) = £o € U the time derivative V = DVf satisfies 
V(x(t,29)) < 0, then w(x) C QNV~1(0). 


Proof. By boundedness of the orbit, w(x) is nonempty and for p € w(x) 
there exists a tg — oo such that a(t,) > p. Since V(x(t,)) < 0 the 
sequence {V(t;,,)} is nonincreasing. Since z(t, £o) is bounded, V (a(tz, £o)) 
is bounded, so that there exists a c € R such that V(t) > c. Hence 
w(x) C V~+(c). Since w(xo) is invariant, w(xo) C Q, and for any y € w(z0) 
we have V (z(t, y)) = c and differentiating gives V (x(t, y)) = 0 for all t, and 
hence V(y) = 0 for all y € w(x). Hence w(x) c QNV~1(0). 


Example 3. 
bony —a(e? +42) 
ý =x +y- yle? +y’). 
Taking U = R?, V (x,y) = y £2 + y2 we get 


dV b 
T =v- vaf 


<0 for |(æ,y)|>1 
> 0 [(2,y)| < 1. 


Thus V~1(0) = {(0,0)} US (S is the unit circle). Q = R? and applying 
LaSalle’s invariance principle we get w((z,y)) C {(0,0)} US. But the 
omega limit set is also connected, so that it must be either {(0,0)} or 
(by invariance) all of S. Since (0,0) is unstable, we must have w(xo) = S 
when xo Æ 0 and w((0,0)) = {(0,0)}. 
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y=ar—(ye+d)y (a, 8,6 > 0) (6) 


This system has a unique steady state (0,0), and one can show that U = 
R&, is forward invariant. Adding (5) and (6) we obtain 


d 
qt) = —dy < 0 on R39. 


Take V(x,y) = £ +y. Then V—1(0) = {(s,0) : s € R}. By LaSalle’s 


invariance principle, 


w((x,y)) E {(s,0):s € Roo}, (x,y) € RSo. 
But w((x,y)) must be connected and invariant, and the only invariant sub- 
sets of T = {(s,0) : s € Ryo} for the flow of (5) and (6) are the origin and 


T itself. But, by (5), for s > 0, ye, (s,0) —> (0,0) for any sequence tg — oo, 
so w((x, y)) = {(0,0)} V(z,y) € R&o. 


Theorem 5 (Goh”). Suppose that the Lotka-Volterra system ti = 
Lif (x) = ari, + Dey aij£j), i = 1,...,n has a unique interior steady 
state x* = —A~'r € RZ). Then this steady state is globally attracting on 
R&, if there exists a diagonal matriz D > 0 such that AD+DA" is negative 
definite. 


Proof. Let V : R$, + R>o be defined by 


n 


V(x) = $ a; (a: — af — 27 log(z:/2})), 


i=l 


where a; € Ryo are to be found. Then we compute 
V=W -f =>) olti- ci) file) = So aile 27) 4 X ailes- 24) 
i=l i=l j=l 
This can be rewritten as 
; 1 
V =(@-2")' AT D(z — z*) = zC — x*)T (DA + AT D) (z — z*), 


where D = diag(a1,..., an). When DA + ATD is negative definite, V < 0 
and V~1!(0) = {a*}. V is convex (as the sum of convex functions) and 
has a unique minimum at x = 2*. Hence by Theorem 3 x* is globally 
asymptotically stable on R&o. 
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(See reference 16 for an improvement of this result to cater for boundary 
steady states.) 


Example 5. Consider the two species Lotka-Volterra system 


t = x(a + bx + cy) 


ý = y(d+ ex + fy). (7) 


Suppose that (7) has a unique interior steady state, say (2*,y*) € R&o. 
Thus bf — ce #0. We use Theorem 5. Let A > 0 and 


10 2b c+ re 
D= M=DA+A™D= f 
o ar a Gar 2f ) 


Then for diagonal stability we need M to negative definite, which it is if 
and only if its trace is negative and its determinant is positive: 


(i) Af +b <0, (ii) 4bfA > (e+ àe)’. 


Since we seek À > 0, to satisfy (ii) we require fb > 0, which then implies 
f,b <0 by (i). Next for (ii) we need 


4bfX— (c+ àe)? = (4bf — 2ce)\ — c? — eA? > 0 


for some À > 0. The quadratic ¢(A) = (4bf — 2ce)\ — c? — e? A? is negative 
for \ = 0 and large |X|, and so is positive for some Aà > 0 only if 4bf — 2ec = 
2det A+2bf > 0 and (4bf — 2ce)? > 4e?c? which simplifies to det A > 0. 


To conclude, we have shown 
Theorem 6 (Goh®). Suppose the system 


t = x(a + bx + cy) 


ý = y(d+ ex + fy). (8) 


has a unique interior steady state (2*,y*) € R29. Then (x*,y*) globally 
attracts all points in R2, if f <0,b <0 and det A > 0. 


Problem 3. Is the converse of Theorem 6 true? 


5. Conservative Lotka-Volterra Systems 


Definition 13 (Conservative Lotka-Volterra). We will say that (1) is 
conservative if there exists a diagonal matrix D > 0 such that AD is skew- 
symmetric. 
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Notice that if B is skew-symmetric then b;; = —b,, for all 7,7. In particular 
bii = —b; so that bii = 0, i.e. the diagonal elements of a skew-symmetric 


matrix are all zero. 


Problem 4. Consider the two-species Lotka- Volterra system 


PN es 
Nda | 

1 dP 

S S NT 
Pa | 


Change to new coordinates x = log N,y = log P and show that H(x,y) = 
dx + ay — e” — e” is constant along a trajectory (x(t), y(t)). Show also that 


A change of coordinates y; = x;/d; (di < 0) transforms (1) into 
n 
Ui = alta + D> dyaiyys), 
j=l 
so that we obtain another Lotka-Volterra system with interaction matrix 
AD. The Lotka-Volterra systems with interaction matrices AD for D > 0 
diagonal have topologically equivalent dynamics. 


Lemma 3. If A is an n x n skew-symmetric matrix then det A = 
(—1)" det A. Hence when n is odd, A is singular. 


Proof. det A = det AT = det(—A) = (—1)" det A. 


Now suppose that A is skew-symmetric. We will show that certain 
Lotka-Volterra systems can be written in Hamiltonian form. But before 
doing so, we recall the definition of a Hamiltonian system on R” (see, for 
example, reference 12). Let C% denote the space of smooth functions 
R” > R. 


Definition 14 (Hamiltonian system on R”). A Hamiltonian system 
(on R”) is a pair (H,{-,-}) where H : R” — R is a smooth function, 
called the Hamiltonian, and {-,-} : C% x C® — C™® is a Poisson bracket; 
that is a bilinear skew-symmetric map {-,-} : C° xC% — C@ that satisfies 
the following relations for all f,g,h € C™ 


(1) {f,gh} ={f,g}h+g{f, h} [Leibnitz rule] ; 
(2) {feta Rh} + {9, {h, f}} + th, {f,9}} =0 [Jacobi Identity] . 
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For example, when n = 2 the bracket {-,-} : C% x C® —> R given by 


Of Og _ Of Og 
19} = ag Bp dp Og 

q OP Pp oq 
defines a Poisson bracket. 

For each g € C™%, the bracket defines a Hamiltonian vector field X, 
on R” via {f,g} = X,(f). In the previous example X} = ee — ap 
Hamilton’s equations are then given by č; = Xp(z;) for i = 1,...,n. 
In particular, H = {H, H} = 0 gives the constancy of the Hamiltonian 
function along an orbit. In addition to conserved functions conserved on 
orbits, there may also be functions C such that {C, f} = 0 for all functions 
f e Cs. That is: C is constant along all flows generated by the Hamilto- 
nian vector fields Xf as f ranges through C'°. Such functions C are known 
as Casimirs. 

To establish that a Lotka-Volterra system is Hamiltonian, we thus have 
to identify both a Poisson bracket and a Hamiltonian function. 

Before turning to a Hamiltonian description of (1) we note that there’s 


a graphical way to test whether a Lotka-Volterra system is conservative: 


Proposition 1 (Volterra’®). The Lotka-Volterra system či = xj(ri + 
(Az);) is conservative if and only if ai = 0 and aij # O > dijaji < 
0, and for every sequence %1,12,...,%5 We have Qijisđizizs''* Qis = 
(—1)°@i,4, 1° ** Qizi Airis: 


That is we have a graphical condition that there exists a diagonal matrix 
D > 0 such that AD is skew-symmetric (AD + DAT = 0; compare with 
Theorem 5). One creates a signed digraph with nodes labelled 1 to n where 
n is the number of species and puts on each directed edge linking nodes i 
to j the number a;;. The condition to check is then that for each cycle in 
the digraph of length s, the product of the edge numbers in one direction 
is (—1)* times the product in the opposite direction. 


6. Volterra’s construction® 18 


We start with the skew-symmetric system 
n 
ti = xilri + Da aijX;), Qij = —ūji. (9) 
j=1 
Volterra introduced new coordinates which he called quantity of life: 


Qi= f alsyas (i=1,...,n). 
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Thus Q; = a; and (9) becomes the second order system 


Q: = Qi(ri+ X aQ). (10) 


j=1 


Then he introduces H (Q, Q) = 37"_, (r:Qi — Qi) so that 


“i =X (riQi-Qi) = X (rQi-Qilr+¢ aQ) =- X yQ = 0 
j j=l ij=l 


i=l i=l 


using skew-symmetry of A = ((a;;)). Dual variables P; are then defined via 


: T , 
P, = log Qi - 3 9 Q; m=) 
j=l 


In terms of these new coordinates, we get the transformed h(Q,P) = 
H(Q,Q) where 


n 1 n 
h(Q, P) = X. | r:Qi— exp(P; + 5 S aQ) 
i=1 j=l 
Now we can check that 
dQi lL ðh 
dt OB wg) OP.’ 


and 


dP, d , IŠ 
= g 42823 N ayQ}; 
j=1 
O; 1 n ; n ; 1 n A 
=-= 32 Q; =r + yQ- 52 aQ; 
Qi j=1 j=1 j=1 
iL ig 
= 7p 3 5 Qij exp(P; + 3 5 ajkQk). 
j=1 k=1 


On the other hand 


Oh R Aki 1 7 7 Qik 1 “ 
m > 5 exp Pat 32 sie Dae a a ’ 


using Qik = —ax;. This gives P; = A as required. 


April 28, 2016 12:43 ws-rv9x6 Book Title LTCClIntegratedBaigent 
page 13 


Lotka- Volterra Dynamics 13 


Hence we have shown that the system (9) is canonically Hamiltonian in 
the new coordinates P,Q with Hamiltonian function 


n 


MP, Q) = 5 riQ; — exp(P; + Sag) 7 
j=l 


i=l 


and the standard Poisson bracket 


Za 
(hah = So Gh oe Se OF 
i=1 


OP,0Q; Q; OP, 


7. An alternative Hamiltonian formulation 


In the previous formulation, we doubled the number of variables in order to 
find a Hamiltonian structure. Here we keep the same number of variables 
as the original Lotka-Volterra system. 

Suppose that Av +r = 0 has a solution z* € R” (here A is skew- 
symmetric). Introduce new variables y; = log x;: 


n n 
Ys = (ri +) aij exp Yj) = X aij (exp yj = v5). 
g=l 


j=l 
Now define 
H(y) =) (expyi — zžyi), 
w=1 
so that 
bs ee 5. ate ". OH 
n= Dale) = aug a) 
j=l ji OMS 
dH QOH. 
de T 2s ay” 
jar 79 
= ina 
i=1 j=1 Oyi OY; 
1S D iy von 
— 9 Qij T Aji) a Du. 
E i=1 j=1 Oyi OY; 
=0 


using skew-symmetry of A = ((a;;)). To complete the Hamiltonian formu- 
lation we check that 


{fg} = Vf- AVg (12) 
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provides a suitable Poisson bracket. 
Problem 5. Show that (12) defines a Poisson bracket. 


Notice that the x* need not lie in the first quadrant. In an odd dimensional 
Lotka-Volterra system with skew-symmetric interaction matrix A, we have 
det A = 0 and it is possible that Ax +r = 0 has no solutions. Indeed, if 
A is singular, then there is a v Æ 0 in ker A such that vT A = (A7v)? = 
—(Av)’ = 0. Thus for a solution to exist we must have vfr = 0 for all 
v € ker A, i.e. r € (ker A)t+. 


Example 6. Consider the Lotka-Volterra system for 3 interacting species: 
tı = zı(rı + wit — W203) 
t2 = z2(r2 — wızı + W323) (13) 
t3 = x3(r3 + W201 — W322) 
where w1,w2,w3 > 0 and each r; > 0. Here species 3 is prey to species 
2. Species 2 consumes species 3, but is consumed by species 1. Species 1 
consumes species 2 but it is consumed by species 3. (So we have a cycle of 
interactions.) It is easy to see that the interaction matrix 
0 Wy, —W2 
A= Wy 0 W3 
wz —w3 0 
is skew-symmetric. Since A is 3 x 3 we already know that A is singular. 
Thus if q is a solution to Aq+r = 0 then so too is q+ k for any k € ker A = 
{a(w3,w2,wW1): a E R}. One finds that Aq +r = 0 has no solutions (in R3) 
unless 
vTr = wer, +wir3 + wore = 0 (14) 
(v = (w3,W2,W1)) and in this case q = (2,—2,0) + av fora ER. 


Wy ? 
Thus let us now assume that (14) holds. For the Hamiltonian we may 


take 
A(z) = 21 + £2 + 23 = log zı + a log £2. 
Wy Wy 
We find that 
Wy Wy 


by virtue of (14). A suitable Poisson bracket is thus 


Ff. g} <urerey (OF 29 _ 29 AF 
i ae we Ox 0x2 Ox 0x2 


Sse Of Og — Og OF | Of Og ög Of 
aes Ox, 0x3 Or, 0x3 phere Ox 0x3 Ox 0x3 i 
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Since A is singular, there are Casimir functions C; that is C satisfying 
{C, g} = 0 for all g, proportional to 
C(x) = w3 log #1 + we log z2 + w1 log z3. 
(or we could take C(x) = xf%a5?a3!). We find that 
C= 73W1 + r2w2 + riw = 0, 

again using (14). The dynamics lies on the intersection of the surfaces 
H (x) = H(«(0)) and C(x) = C(x(0)) in the first quadrant. 

Let us change coordinates, setting X = loga1, Y = logx and Z = 


log x3. Then we have on a solution 


rE eee ee re a 
Wy Wy 


w3X + wY + w1Z = B, 
where A, B are constants. Hence we may plot 
T2 


Z = log (4 eae ae a Y) (15) 
WL Wy 


= B-w3X — wY 
= m . 


Z (16) 


The first surface is concave where the logarithm is defined. Searching for 
periodic orbits then becomes the study of how the surface (15) intersects 
the plane (16). An example a periodic orbit is shown in Figure 1. 


8. Cooperative Lotka-Volterra Systems 


We will consider the general Lotka-Volterra system 
ti = F(x) = xi(ri + X agjay), (@=1,...,7). (17) 
j=l 


except that we will constrain ourselves to the case that a;; > 0 when i Æ J, 
i.e. the off-diagonal elements of the interaction matrix are nonnegative. 
Notice that in this case 
OF; 
Ox; 
since for i # j we have aj; > 0 and we have x € RZ. Since the first 
quadrant is invariant the Jacobian has nonnegative off-diagonal elements. 


= ijti > 0, iF fj, 


Definition 15 (Cooperative matrix). We will say that any real n x n 
matrix with nonnegative off-diagonal elements is cooperative. 
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Fig. 1. A periodic solution to the three species model (13). There exists a continuum 
of periodic orbits around the interior steady state. 


Some notation 


In what follows we will use the following notation for ordering vectors x € 
R”: For each x,y € R” 


e r< y &ri< yi for alli=1,...,n; 
e r<y $ ri< yi for alli=1,...,n, but £k Æ Yk for some k. 
oe r&y ri < yi for alli=1,...,n. 


(Similarly for >, >, >.) We say that (R”,<) is an ordered vector space. 


The following Perron-Frobenius theorem is fundamental in the study of 
coorperative or competitive systems (see, for example, reference 3). We 
recall that the spectral radius of A, written p(A), is the modulus of an 
eigenvalue of A of largest modulus, and a matrix A is irreducible if it is 
not similar via a permutation to a block upper triangular matrix (that has 
more than one block of positive size). 
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Theorem 7 (Perron-Frobenius). If A is an xn real matrix with non- 
negative entries. Then 


e p(A) is an eigenvalue of A. 
e A has left and right eigenvectors u > 0 and v > 0 associated with p(A) 
(i.e. uA = p(A)u and Av = p(A)v). 


If A is also irreducible then we have p(A) > |p| for any eigenvalue u 4 p(A), 
and p(A) is simple and v > 0 and u > 0 in the above statements. 


(Note: the inequalities in Theorem 7 use the vector ordering defined above). 


Definition 16. A matrix A is negatively (row) diagonally dominant if there 
exists a d > 0 such that andi + $ 54; aijld; < 0 for all i =1,...,n. 


When A is a cooperative matrix this becomes Ad < 0. 


Lemma 4. Let A be a cooperative matrix. Then A is stable if and only if 
it is negatively diagonally dominant. 


Proof. First suppose that A is negatively diagonally dominant: There 
exists a d >> 0 such that Ad < 0. Note that we must have all a;; < 0 
since the off-diagonal elements are nonnegative and d > 0. Let A be an 
eigenvalue of A with right eigenvector x. Let y; = xi/di for i = 1,...,n 
and |Ym| = max; |y;| > 0. Then Adiyi = 2a aijdjyj and 


Adm = dm Grim + y Dyamy t 


jem É 
Therefore 
n Yj n 
[Adm = dm@mm| < 5 djQmj u < DD djQmj < —dmamm 
m 


jxm jxm 
by hypothesis. Hence |À — amm| < —amm and A must lie in the open disc 
in the Argand plane whose boundary passes through zero and whose centre 
is at the negative number amm. Thus all eigenvalues \ have negative real 
part, so A is stable. 

Conversely, suppose that A is stable and has nonnegative off-diagonal 
elements. For c > 0 sufficiently large B = A + cI is a nonnegative matrix 
and so by the Perron-Frobenius theorem there is a A = p(B) > 0 anda 
v > 0 such that Bv = Av = p(B)v. But then Av = (p(B) — c)v so that, 
since A is stable, p(B) < c. Since p(B) < c the following series converges 


ee G ‘BY aE +: 
Cc Cc 
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and thus all elements of A`! are non-positive. Now set d = 
—A-1(1,...,1)7. Then d > 0 (no row of A can be zero, since it is nonsin- 


gular) and Ad = —(1,...,1)7 <0. 


As a corollary we have: 


Corollary 1. If A is cooperative and r >> 0 then Ax +r = 0 has a unique 
interior solution x € R& if and only if A is stable. 


Problem 6. For the system (8) when c > 0, e > 0 and det A > 0, find the 
condition for a unique interior steady state. 


We also have the following (see, for example, Theorem 15.1.1 in reference 
10): 


Theorem 8 (Global convergence for cooperative Lotka-Volterra). 
When A cooperative and stable the system (17) has a unique interior steady 
state that attracts R&o. 


Proof. By corollary 1 A is negatively diagonally dominant by lemma 4, i.e 
there exists a d >> 0 such that a;;d; + Bia aijd; < 0. Define 


[Ek — Thl 
V(x) = ——, 
(x) max ak 
Then V(x) > 0 with equality if and only if x = x*. Now, consider a time 
interval during which max; Te = ee, Then 
V = —tisgn(zi — 27) 
di 
Ti È * * 
= a ij (Xj — UF) +X a(z; — %5) p sgn(z; — £} ) 
ý j+i 
Ti ok * 
< ET anle: — efl + X ailt; — 27 
di i#i 
ety d d 
S (x) T S ot) 


< 0 for all x € RZo, with equality if and only if x = x*. 


Hence by Theorem 3, x(t) > x* as t > œ. 
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9. Competitive Lotka-Volterra Systems 


Now we consider the Lotka-Volterra system 
n 
j=l 


under the special conditions that aj; > 0 for all 1 < i,j < n (caution: 
notice the change of sign in (18)). This means that each species competes 
with all other species including itself. If some r; < 0 then it is clear that 
x(t) > 0 as t > œ since R&, is invariant and 


n 
X i 2 
Ti = zilri = Ai; 25) < hiiti < 0, 
j=1 


with equality if and only if x; = 0. We will therefore also assume r; > 0 for 
each i = 1,...,n. This means that in the absence of any competitors the 
species i will evolve according to t; = x;(r; — aixi) and hence will either 
remain at zero or stabilise at its carrying capacity K; = r;/a;; > 0. It also 
means that the origin is an unstable node. 


Lemma 5. Since aij > 0 and r; > 0, all orbits of (18) are bounded. 


Proof. R&_ is invariant and 


n 
r; 
r 2 ; i 
Zi = Tifi — Ti y Qij£j LX; — Quti = zilri = Qizi) < 0 if x; > To 
Ñ a 
j=1 


so that the ith species is bounded for each i = 1,...,n. 


10. Smale’s Construction 


In the 1970’s many thought that for a finite habitat that is home to a 
number of species that compete with each other and the other species, the 
long term outcome is “simple” dynamics, e.g. convergence to a steady state 
or a periodic orbit. But this is not the case, as Stephen Smale showed in 
1976.14 Consider a more general model of total competition: 


where M; is smooth and we will suppose that 


S1 For all pairs i,j we have h < 0 when z; > 0 (totally competitive). 
S2 There is a constant K such that for each i, M;(x) < 0 if |a| > K. 
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Condition S1 means that 
OL; OM; se sh 
“i = <0 alli,j if z; >O. (20) 


dx; | Ox; 
Thus the Jacobian has negative elements in R&,. In other words competi- 
tion for resources. The second condition says that there are finite resources 
and that the populations can not grow indefinitely. 

Smale showed that examples of systems satisfying (19) and the condi- 
tions S1, S2 whose long term dynamics lie on a simplex and obey ¢ = h(x) 
on the simplex, where h is any smooth vector field of our choice! Thus the 
simplex is an attractor upon which arbitrary dynamics can be specified, 
even chaos. 

We follow the presentation of reference 9. Let Ay = {x € R&q: lell = 
1} be the standard probability simplex with tangent space Ap = {x € R” : 
X; vi = O}. Let ho : Ay —> Ao be a smooth vector field on A; whose 
components can be written as hoi(x) = xgi(x) and h : R&, > Ao any 
smooth map which agrees with ho on Aj. 

Now let 8 : R > R be any smooth function which is 1 in a neighbour- 
hood of 1 and A(t) = 0 ift < 4 ort > 3. For e > 0 define M; on R&y 
by 

M,(x) =1— [lalla + (llel )gi(), 1 <i <n. 
We may check: for each i, j, 
OM; 
Ox; 
for small enough e since @ has compact support. 
Now as before, R&, is invariant, and £ |æ] = Xi t = |z- llel) 


(the logistic equation!). Thus A, is forward invariant and any point in 
RZ, \ {0} is attracted to A;. On A; we have 


M,(2) =1— lalla + 8 (ell) 9i(@) = ege(2), 


so that the dynamics on the attractor is t; = xcgi(x) = eh;(x) for i = 
1,...,n, with h arbitrary. 

Hence we should be warned that the long term dynamics of bounded 
competitive systems in dimensions higher than two can be very complex 
(although one can show [see the next section on the carrying simplex] that 
when n = 3 the long-term dynamics must lie on a set of dimension at 
most 2, and this severely restricts the possibilities. However, much more is 
possible when n > 4.) 


Ogi 
Ox; 


= —1 + ep (|læll1)gi + B(l2|11) <0, 
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11. Carrying Simplices 


A bounded totally competitive system with the origin unstable has a unique 
invariant manifold that attracts the first orthant minus the origin. We will 
give an example* of such a system where the invariant manifold can be 
explicitly found - it is the probability simplex in R&, - and all orbits save 
the origin are attracted to it. Moreover, (for that example) the dynamics 
on the simplex is canonically Hamiltonian and all orbits are periodic. 

We consider again the system 


Xi = xiıMı(zx), (i a Ls eN n), 


where M; is smooth and we will suppose that 


S1 For all pairs i,j we have h <0. 
% 


S2 There is a constant K such that for each i, M;(x) < 0 if |x| > K. 
$3 M,(0) > 0. 


Condition 53 makes the origin 0 a repelling steady state. Since orbits are 
bounded, the basin of repulsion of 0 in R&, is bounded. The boundary of 
the basin of repulsion is called the Carrying Simplex and is denoted by X. 
One can think of © as being the boundary of the set of points whose a limit 
is the origin. 

All steady states and all w limit sets lie in © and we have from Hirsch® 


Theorem 9 (The Carrying Simplex). Given (19) every trajectory in 

Zo \ {0} is asymptotic to one in X, and X is a Lipschitz submanifold, 
everywhere transverse to all strictly positive directions, and homeomorphic 
to the probability simplex. 


Thus totally competitive n—dimensional Lotka-Volterra systems (as above) 
eventually evolve like n — 1 dimensional systems. Thus nothing very exotic 
can happen for n < 4. In Figure 3 we display 3 examples of the carrying 
simplex for totally competitive Lotka-Volterra systems. 

The following example has the advantage that the carrying simplex can 
be found explicitly, and it is easy to see that all points save the origin are 
attracted to it. 


Example 7. We consider the illustrative example of an eventually periodic 


4A second example, since in Smale’s example the unit simplex is also a carrying simplex. 
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2 


contains the 
steady states 


Fig. 2. The Carrying Simplex attracts all orbits except the origin and contains any w 
limit set and in particular all steady states except the origin. 
competitive system! 
t = a(1—x2- ay — Bz) 
ý = y(1 — Br — y — az) 
ż = z(1 — az — By —z) 
where a, 8 > 0 and a + 8 = 2. Let 


V(£,Y, Z) = ryz. 


Then 
GV = ay: (24242) 
=V ((1 -z — ay — bz) + (1 — px — y — az) + (1 — az — by — z)) 
=V(3— (x +y +z) -— (a+b) +y+z2)) 
= 3V(1—(x+y+z)) since a +8 =2. 
Moreover 


L e+y+) = (@ y +2z)—= z — y — z’ — (a+ B)(ry + grz 4+ yz) 
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Examples of the carrying simplex for competitive the 3 dimensional Lotka- 


Fig. 3. 
Volterra equations. From left to right the carrying simplex is (i) convex, (ii) concave 


and (iii) saddle-like.! 

Thus if (0, yo, zo) € R3 \ (0,0, 0) we have x(t) + y(t) + z(t) 3 1 as t > o. 
That is, all orbits eventually end up on the simplex A;. Thus the carrying 

simplex © in this example is just the simplex A,;. On A, we have 


dt 
that is V = const on A;. What is the dynamics actually on the carrying 


simplex? We may eliminate z since z = 1 — x — y on the carrying simplex. 


dV 
3V(1—(x+y+z))=0, 


This gives 
z-ay- 60 -e-y) = CZP- 2-2) 


a(l-2-y))=—% 


t= x(1 


ý = y(1 - px- y 
where a + 8 = 2. Notice that div («,y) = 0 and that we have a canonical 


Hamiltonian system with Hamiltonian function 
(a— 6) 
Seay. 


A(x,y) = 2 


On the open triangle T = {(z,y) € Rọ: 0 < x+y < 1} we obtain 
closed contours, i.e. the solutions are periodic. (This is the projection of 
the dynamics on © onto the zy—plane.) Figure 7 shows the periodic orbits 
on the invariant plane © = {x € R3% : 21 + £2 +23 = 1} as part of the 


3-dimnesional phase portrait. 
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Fig. 4. Periodic orbits in a model of May and Leonard.!! Note the carrying simplex is 
the usual simplex in Ro and it clearly attracts all orbits apart from the origin. 


Problem 7. Consider the planar competitive system?” 


é y 
= 1 m TE 
z= x( x 5) 


ý = y(1 — 3x — y). 


Show that this system has a carrying simplex © which is the graph of a 
quadratic function q that satisfies q(0) = 1 and q(1) = 0. Sketch the phase 
plane. 


12. Further reading 


The book by Takeuchi (reference 16) provides a comprehensive study of 
Lotka-Volterra equations. Hofbauer and Sigmund!’ contains much on 
Lotka-Volterra systems, plus also their close cousins the Replicator equa- 
tions from evolutionary game theory. The cooperative and competitive 
Lotka-Volterra models discussed here are a small subset of monotone dy- 
namical systems which are covered in reference 9 by Hirsch and Smith. See 
also the monograph on monotone dynamical systems by Smith.!° In 2003 
a new geometrical approach to the study of Lotka-Volterra systems was 
initiated by M. L. Zeeman and E. C. Zeeman.!® The carrying simplex of 
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competition models originated in a 1988 paper by Hirsch. The geome- 
try of carrying simplices has been studied by Tineo!” and more recently 
by Baigent,' where in the latter paper links are made to the geometrical 


approach of Zeeman and Zeeman.!* 


13. Sketch solutions to problems 


Problem 1: 
&=a(1—x2-Yy) (21) 


y=y(-l—-y+2). (22) 
From (21) we see that for any initial condition, since y(t) > 0, for t large 
enough 0 < a(t) < 1. But then, from (22) eventually y < 0 and so we 
must have y(t) > 0 as t > oo. It is thus clear that all orbits are bounded 
and each w((xo,yo)) is nonempty. If (p1,p2) E w((xo, yo)), then for some 
tk > co, po = limg+co y(t) = 0. Thus w((xo0,yo)) C {(s,0) : s € [0, 1]}. 
If (Py(t),0) is the forward orbit with (P,(0),0) = (91,0) E€ w((2o, yo)) 
and pı > 0 then by invariance (P\(t),0) € w((xo,yo)) for all t € R. 
But P,(t) satisfies È, = P,(1 — Pı) and so py = liMmk>o Piltk) = 1. 
Hence for any (£o,yo) € Ro \ {(0,0)} we have w((xo,yo)) = {(1,0)}, 
and w((0,0)) = {(0,0)}. 
Problem 2 
The system t = x(1 — x) has invariant sets {0}, {1}, [0,1], {0,1}, Rso, 
R>o \ {1}. The forward invariant sets are [s1, s2] for any 0 < sı < 1 < s2. 


B(1) = Ryo. 

Problem 3 

No, the converse does not hold. Take the system t = x(1— § — $) and 
y = y(-1+a+ 4). Then there is a steady state ($, $) which globally 
attracts R29, with det A = 4 > 0 and b = —4 < 0 but f = $ > 0. 
Problem 4 

Set « = log N,y = log P. Then « = a — be”, ġ = ce” — d. Set H(z, y) = 
e Then oH — = d— ce? = —y and on =a—be’ = t. Then 
H= ait Hi +4 ty = = 0, so that H is constant along ae 

Problem 5 


{f,g} = Vf - AVg, where AT = —A. Then using the summation con- 
vention {g,h} = aij oe oe . Thus {f,{g,h}} = aij sL F (an n 2) = 


Tað 3g Oh, Og _&Æh \_ Of dh g , Of Og Əh 
Qij Qik Hi Ox, | Ox; Ox,02; ) — Tijk \ Ix; Day Orð; | Dx; Ox; Ot,02; | 


: _ Og Of Əh Og Oh Əf 
By cycling terms {9, {h, fh} = ajan ($2 xk OL1ðTj | Ox; Ox, OL,OL; 
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cae dh 0g Of adh Of dg ; f 
and {h,{f,g}} = Qija (2 Oa, Ones Ot: Ozi Dero, J? Consider sec- 


ond derivatives of h the sum {f, tor) g 19, {h, f}} + {h, Vy. Dy this 
Of Og Əh 


| = —_ 
gives Qik Jr; Bx, Ox,0r; | Qijalk af A aeos = GjUk dz, ne LA 
dg Of h _ ah_ (dg Af _ af dg 
ij Uk Jr; Oa; Ox,0a; — VINK Fz, 0a, (#2 Tx; — Der De, | This value of this 


last expression is not changed by interchanging labels 7,1, whereas the 
bracketed term changes sign, and so the expression must be zero. 
Problem 6 

The conditions are det A = bf — ec > 0 and dc > fa, ea > db. Since ec > 0 
we must have bf > 0, and hence either (i) b, f > 0 or (ii) b, f,< 0. If 
b, f > 0 then 7d >a and ţa > d. From these it is clear that both cannot 
be satisfied unless a,d > 0. If b, f > 0 and a,d > 0 then the condition is 
8 <ISF (Of course, even if an interior steady state exists, the cooper- 
ative system (e, c > 0) may have unbounded orbits). 

Problem 7 

An invariant curve connecting (1,0) and (0, 1) is a solution y : [0,1] + R>o 
of y' (x) = ee that satisfies the boundary conditions y(0) = 1 and 
y(1) = 0. Let the quadratic be y(x) = 1+ ax + bz?. Then this sat- 
isfies the boundary condition y(0) = 1. To satisfy y(1) = 0 we need 
0 = 1 +a +b so that we need a = —1 — b, and y takes the form y(x) = 
1—(1+b)x+ bx? = (1—2)(1—bz). For this function y' (£) = —1 — b+ 2bz, 
so that y’(0) = —1 — b. On the other hand y'(x) = YU=32—) so by 


a(l-a—¥)? 
L’H6pital’s rule y/(0) = ole S © Now set y(0) = 1 


to obtain y’(0) = —2, which implies ie 1, and thus that the curve is 

y(x) = (1 — x)?. It is now simple to check that y(x) = (1 — x)? satisfies 

y' (x)= ea and so the graph of y is an invariant manifold. See Figure 
2 


5 for the phase portrait. 
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